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The effect of nonmagnetic disorder on the pairing amplitude is studied by means of a 
perturbation method. Using an extended one band Hubbard model with an intersitc 
attraction we analyze various solutions of p-wave pairing symmetry and discuss their 
instability calculating fluctuations of order parameter. The model is applied to describe 
the effect of disorder in Sr2Ru04. The results shows that the real solution with line 
nodes can be favoured by disorder. 
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1. Introduction 

The recent discoveries of superconductivity in Sr 2 RuO 4 0§ and UPt JB, apparently 
p-wave, raised the questions of their unconventional properties. Particularly the 
question of svstom r e s p onse on nonmagnetic disorder attracted interest of many 
researchersuu-B0l§Bll3liJ . .It-is clear that this effect can distinguish the uncon- 
ventional superconductorai3'L2l. For conventional superconductors (with is otropic 
s-wave paring and a large coherence length) it is known, due to Andersonli^tJO 
that the non-magnetic disorder have negligible effect on superconducting gap and 
the critical temperature Tq- 

The fluctuations oLthe pairing amplitude are scaled as the inverse of coherence 
length square (~ £ _2 )li3. In an analogy to that, for anisotropic supercojuiuctors, 
disorder induces the effect of spatial fluctuations in the amplitude | A(zj)|EJ3c3 which 
can destroy the superconducting phase. 

In this paper we will consider the effect of disorder on various possible p-wave 
solutions. It is connected with a dilemmaEl appeared recently in context of super- 
conductivity in Sr2RuC>4. Namely, the question if superconducting order parameter 
posses or not the line nodes in the gap was asked. As among possible p-wave paired 
solutions there are those with line nodes and without them, therefore it is important 
to check if disorder can favour any of possible p-wave solutions. 

For simplicity we are going to use the one band modclESS. Such choice can 
be justified knowing that one of the promising model of p-wave superconductibity 
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in Sr2Ru04 is based on the important assumption that 7 band can be dominant 
as far as the appearance of superconductivity is concerned while the electrons in 
other two passive bands (a and 0) are paired because of proximity effectE3. 

The paper is organized as follows. After short introduction in the present sec- 
tion we present the model in the clean limit and introduce the Hamiltonian and 
approximations used in paper (Sec. 2). Here we discuss the results of paring pa- 
rameter obtained for triplet superconductor. Section 3 is devoted to a disordered 
p-wave superconductor. We include diagonal nonmagnetic disorder and analyze 
spatial fluctuations of a pairing potential. Particularly we calculate standard devi- 
ations random site energies < fe; > and for pairing potentials for neighbour sites 
(ij) < 5Aij >. Their ratio =j£_5Ajj > / < 5ei > can be regarded a criterion 
of pairing potential fluctuationsMLia Finally we investigate the stability of various 
solutions in presence of disorder. Section 5 contains conclusions and remarks. 



2. The model. 

We start a single band, extended, Hubbard model with intersite attraction. It is 
defined by the following HamiltoniarcJ: 

H = ^{e l 8 t j + t i:i )cfc ja + -^-rv i<T n j(Tl - n^c+da , (1) 

ija ijacr' ia 

where cf a and ci a are usual, fermionic operators which create and annihilate, re- 
spectively, an electron with spin a at the lattice site labeled by i, tij is a electron 
hopping integral, is the operator of particle number of spin a at site i, e< is 
the site energy, varying from site to site in random fashion, at the site i with mean 
value £0 — < £j >— and Wij is the interaction potential of two electrons with 
opposite spins on neighbour sites /1 denotes the chemical potential. 

Hartree-Fock-Gorkov equation of motion for matrix Greens 4x4 functions 
G{i,j;uj) yields: 

E( [(icon - £i + fi)5 a + in] 1 A a \ (Gu(l, i;iu n ) Gi 2 (l,j;iu> n )\ _ r 1( .„s 
. ^A+ [{iuj n + e l -^5 a -t l i]l)\G 2 i{l,j^ n ) G 22 {l,j;iu n )J ~ ij U 

with the spin dependent 4x4 Green function. Each of its components G nm (l,j; ioj n ) 
is defined as: 

Gn ro a,i;^)=fSH! ,i; "ir?H5 ,i; "iy ^=1,2. ( 3) 

\Gim(hj;iuv) G# m {i,3\ih} v ) J 

The order parameter for triplet paring reads: 

/ ATT aTJ- \ 

y = A 1 J T Ji ■ ( 4 ) 
The gap equation for p-wave order parameter: 

A ?/ = U »\Y, ^ nV G^\i,j; iu n ), a, a' =T, 4 , (5) 
P ^ 
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where G"^* is a matrix element ol 4 x 4 functions Green functions. 

The above equation Eq. || have to be completed by the corresponding equations 
for the chemical potential \i that satisfies the following relation: 

n=-Y.^G^{i^n\ (6) 

n 

where n is the number of electrons per unit cell. 




Fig. 1. The Fermi surface (a) and the density of states N(E) (b) for 7 band Sr2Ru04 (t' = 0.45t, 
n = 1.34). 



3. P wave solutions for a clean system at T — OK 



G°(fc;w) 



(7) 



Taking the Fourier transform of the Green function for a clean system (gj = 0) we 
can write its equation of motion: 

(cj-e fc + /i)l A fc 

Here, we assumed that the electron hopping integrals tjj can take a nonzero 
values for the nearest and next nearest neighbour sites only. For a clean system 
can be expressed in fc-space: e k = J^. t i je"*' R « fe as: 

£fc = — 2i(cos k x + cos k y ) — 4t' cos k x cos k y , (8) 

where t represents the nearest neighbour site amplitude of electron hopping, while 
t' corresponds to next nearest neighbour one. 

As can be readily shown the p-wave superconducting order parameter A fc can be 
written as 



A fc = A x sin k x 



A y sin ky, 



(9) 
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where Aj and A. y are complex 2x2 matrix of pairing potentials. Assuming the 
order parameter A fe = ia y a ■ dtkl with d(k) = (0, 0, d z (k))\Q and cP(k) = A k The 
gap equation can be written asE3: 



N ^ 2E q \ 2 J 



where for p-wave pairing we can assume 

W kq = \W\ (2s\n(k x )sai(q x ) + 2sm(k y )sm(q y )) 
and E q denotes quasi-particle energy: 



i 2 - A 2 
q i 



and the equation for chemical potential 
n- 1 



lyls_ Tanh^^ 
N ^ 2E q 



(10) 



(11) 



(12) 



(13) 



The corresponding free energy F for a finite temperature T can be calculated 
from the following formula: 



I k ; , 2k L; T In ( 2 cosh 



|A fc (r)| s 
w 



(14) 



To perform numerical calculations we have fitted our one band system param- 
eters to the realistic 7 band structure of S^RuoJjfflla'El. Fig. 1 presents the 
corresponding Fermi surface (Fig. la) as well as the electron density of states 
N(E) (Fig. lb) of the normal state. Solving Eqs. |l^-|3| we have calculated the 
pairing parameter A x and A y for various possible p-wave superconducting states. 
Namely, depending on relative values A x and A y : dipole (A x ^ and A y = 0), 
real (A x = A y ) and complex one (A^ = iA y ) where found. In Fig. 2a we plotted 
A x for these solutions versus band filling n. The lines are denoted by letters: D, R, 
C; respectively. 

The interaction parameter used in calculations (W = — 0.4605t) were chosen 
to give the same T c — 1.5K for all p-wave solutions, as for the pure Sr2Ru04. 
Examining the free energy F (Eq. |l4|) we have also found the energetic preferences 
of these three solutions. In Fig. 2b the free energy F was plotted as a function of 
|A|. In One can easily see that the complex solution reaches the global minimum 
of free energy F. Additionally, the angular dependences of gaps ((f) S [0, 2n}) have 
been plotted in Figs. 3a, b and c. Note that both dipole and real solutions posses 
line nodes while the complex one has a finite gap in any angle cf> direction (<f) is 
defined in Fig la). For clarity we have plotted corresponding quasi particle density 
of states in Figs. 3b, d and f, respectively. Note that only in case with the complex 
solution it has a finite gap (Fig. 3f). Interestingly the real solution has the largest 
amplitude (Figs. 3c, d). 
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Fig. 2. (a) the pairing parameter A 
of |A| (|A| 



versus band filling n and (b) the free energy F as a function 



A|) for different solutions with dipole (D) real (R) and complex (C) order 

0.4605* and temperature T = OK. The 



parameters, respectively. The intersite attraction W 
arrow in Fig. 2a corresponds to the 7 band filling in Sr2Ru04 n 



1.34. 



4. Fluctuations of pairing potential. 

In this section we investigate the stability of superconducting p-wavestates in 
presence of weak disorder. Here, we apply the same strategy as ir# 6 H 18 H 19 l and we 
treat random site energies £j as perturbations. To proceed we write the Dyson 
equation for a Green function G(i,j;uj) evaluated at a frequency uj: 

G(i,j;u)=G°(i,j;u ) ) + Y,G°(i,l;oj)V l G(l,j;u) , (15) 

1 

where Vi defines the impurity potential matrix: 

V-fa°J. (16) 



-e, 

To the lowest order in Si we get: 

G(i,j;oj)=G°(i,j;oj)+J2G°(i,l;u;)V l G (l,j;u;) (17) 



Following Eqs. (jq,|r7|) we express interesting quantities A^- and Ai/llSta, where I 
and j are nearest neighbours of i, in the lowest order of £j perturbations by means 
of disordered Green function (Eq. |l^) and calculate the mean square deviation of 
the paring parameter as: 

<|JA ii | 2 >=<|A ii | 2 >-|<A ii >| 2 , (18) 

where i and j are nearest neighbour lattice sites. 

The assumption is that random site energies £j in Eqs. (|5| Ji^Ji"8| ) are independent 
variables then averages < Si > ,< £, > being independent of the site index i and 
< £i£j >=< £? > Sij lead to 

< SA^SA*^ >= r, 3 < e\ > (19) 
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Fig. 3. The angular gap (±i?^,(e^, = fi)) dependence and the quasi particle density of states for 
different solutions: dipole (a,b) real (c,d) and complex (e,f) type order parameter; respectively. 
The angle <j> is defined in Fig. la. 



Finally, we calculate the coefficient 1^-00 



1-3 



2N k {Ek+ E k-q) E k E k-q 



he k-q i(R t -Ri)k 



(20) 
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Fig. 4. Fluctuation parameter T as a function of band filling for different solutions: dipole (D) 
real (R) and complex (C) order parameter, respectively. The intersite attraction W = —0.46054 
and the band filling n = 1.34 correspond to 7 band occupation of Sr2Ru04. 

Having solved the gap equations (Eqs. |9r|i"i"|) for the clean system , we calculated 
r (Eq. |2(]) for each of solutions we obtained. They are shown in each of Fig. 4, where 
we plotted T versus band filling n. Analyzing the results, first of all, we have to 
report that for all cases T has small value (of order ~ 10~ 7 ). This implies that 
fluctuations of Ay are small in this system and it is consistent with relatively large 
coherence length in this superconductorll3. Interestingly, a real type of solution has 
the smallest value of T. This could mean that this solution is favorable by disorder. 
Note that the maximum of T (Fig. 4), however close to the van Hove singularity, 
occurs for a lower band filing n than the maximum of the order parameter (Fig. 
2a). 

5. Conclusions and Discussion 

We have analyzed the effect of weak disorder £>n p-wave superconductor in context 
of newly discovered si 1 percon H n ctor Sr&Ru Q Jj 'a. The order parameter structure in 
this compound is still under debateEI'Ei'ESSEa. 

Fitting our system parameters to its 7 band structure we asked about the sta- 
bility of various solution in presence of disorder. Note that the standard Abricosov- 
Gorkov formalismLO applied for site diagonal nonmagnetic impurities could give us 
only partial answer on the disorder effect on superconductivity in Si^RuO,!. It en- 
ables us to predict how strongly the critical temperature T c is suppresedBu but it is 
not useful to analyze the stability of different solutions with the same T c appeared 
in our case. 

Our preliminary results indicate that the complex solution has the global mini- 
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mum of free energy F but it is the real one, with the line nodes, which is favored 
by disorder. That result was obtained in one band model in the lowest order of 
perturbation method and should be investigated further by a more sophisticated 
method which can treat a finite disorder, and by considering more realistic three 
dimensional, three orbital electronic structure of S^RuOj^a. 
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